Paramagnon-induced dispersion anomalies in the cuprates 
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We report the self-energy associated with RPA magnetic susceptibility in the hole-doped 
Bi2Sr2CuOe (Bi2201) and the electron-doped Nd2- : rCe :! ;Cu04 (NCCO) in the overdoped regime 
within the framework of a one-band Hubbard model. Strong weight is found in the magnetic spec- 
trum around (7r,0) at about 360 meV in Bi2201 and 640 meV in NCCO, which yields dispersion 
anomalies in accord with the recently observed 'waterfall' effects in the cuprates. 
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Very recent angle-resolved photoemission (ARPES) ex- 
periments in the cuprates have revealed the presence of 
an intermediate energy scale in the 300-800 meV range 
where spectral peaks disperse and broaden rapidly with 
momentum, giving this anomalous dispersion the appear- 
ance of a 'waterfall' [E 0, S S> H> @] • Similar self-energies 
have also been adduced from optical data0j- This new 
energy scale is to be contrasted with the well-known 
low energy 'kinks' in the 50-70 meV range, which have 
been discussed frequently in the cuprates as arising from 
the bosonic coupling of the electronic system with either 
phononsfH and /or magnetic modes Although low en- 
ergy plasmons 1^, U | are an obvious choice for the new 
boson, analysis indicates that the plasmons lie at too high 
an energy of ~1 eV to constitute a viable candidate [12J. 
Here we demonstrate that paramagnons provide not only 
an explanation of the energy scale but also of the other 
observed characteristics of the waterfall effect in both 
hole and electron doped cuprates. 

For this purpose, we have evaluated the self-energy 
associated with the RPA magnetic susceptibility in the 
hole-doped Bi2Sr2Cu06 (Bi220f ) and the electron-doped 
Nda-zCe^CuCU (NCCO). [13 In order to keep the com- 
putations manageable, the treatment is restricted to the 
overdoped systems where magnetic instabilities are not 
expected to present a complication. Our analysis pro- 
ceeds within the framework of the one-band Hubbard 
Hamiltonian, where the bare band is fit to the tight- 
binding LDA dispersion 14 , HH . We incorporate self- 
consistency by calculating the self energy and suscep- 
tibility using an approximate rcnormalized one-particle 
Green function 



function. The RPA susceptibility is given by 



G = Z/(u-£ k +i6), 
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where = Z(ek — /-*)■ Here, are bare energies and /i 
is the chemical potential, and the renormalization factor 
is Z ~ (1 — dYi 1 / doj)^ 1 < 1. The associated magnetic 
susceptibility is 
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with U denoting the Hubbard parameter. The self-energy 
can be obtained st raig htforwardly from the susceptibility 
via the expression [16} (at T = 0) 
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where S is a positive infinitesimal, fr = /(en) is the Fermi 



Concerning technical details, we note that for the 
generic purposes of this study, all computations in this ar- 
ticle employ a fixed value Z =0.5, which is representative 
of the band dispersions observed experimentally in hole 
as well as electron doped cuprates. [17] Self-consistency 
is then achieved approximately by determining values of 
the chemical potential fi and the Hubbard parameter U 
to keep a fixed doping level and to ensure that the bands 
are indeed renormalized by the average factor Z = 0.5. 
The procedure is relatively simple, but it should capture 
the essential physics of the electron-paramagnon inter- 
action, although our treatment neglects the energy [lijl 
and momentum dependencies of Z. Note also that in 
the overdoped regime considered, the effective U values 
in Bi2201 and NCCO are small enough that the system 
remains paramagnetic and the complications of the anti- 
ferromagnetic instability are circumvented. Specifically, 
the presented results on Bi2201 are for x = 0.27 with 
(jt = -0.43 eV and U = 3.2t, while for NCCO, x = -0.25 
with fi = 0.18 eV and U = it. 

Figure U summarizes the results for Bi2201. We con- 
sider Figs. 1(a) and (b) first, which give the real and 
imaginary parts of the self-energy at several different 
momenta as a function of frequency. The theoreti- 
cal self-energies, which refer to Bi2201, should be com- 
pared directly with the corresponding experimental data 
(gold squares 0), although available experimental points 
for Bi22120 and La 2 _ ;c Sr ;c Cu04 (LSCO)0 are also in- 
cluded for completeness. The agreement between theory 
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and experiment is seen to be quite good for the real part 
of the self-energy in (a) , while theory underestimates the 
imaginary part of the self-energy by a factor of ~ 2. That 
the computed X" is smaller than the experimental one is 
to be generally expected since our calculations do not 
account for scattering effects beyond those of the para- 
magnons. Here, we should keep in mind that there are 
uncertainties inherent in the experimental self-energies 
due to different assumptions invoked by various authors 
concerning the bare dispersions in analyzing the data. In 
particular, Feng et a/.[4| extract the bare dispersion by 
assuming that X' is always positive and goes to zero at 
large energies. Other groups]^, fl^ | compare their results 
to LDA calculations and argue that X' must become neg- 
ative at higher energies. Our computed X' in Fig. 1(a) 
becomes negative over the range 0.35-0.9 eV in certain k- 
directions. Interestingly, various computed colored lines 
in (a) and (b) more or less fall on top of one another, 
indicating that the self-energy is relatively insensitive to 
momentum, especially below the Fermi level, consistent 
with experimental findings H, even though X possesses a 
fairly strong frequency dependence. 

Fig. 1(c) gives further insight into the nature of the 
spectral intensity obtained from the self energy of Eq. [4] 
The spectral intensity shown in the color plot of the fig- 
ure is representative of the ARPES spectrum, matrix 
element effects[2(| notwithstanding. The peak of the 
spectral density function defined by taking momentum 
density cuts (MDCs), shown by yellow dots, follows the 
renormalized dispersion (orange dashed line) up to bind- 
ing energy of about 200 meV. It then disperses to higher 
energies rapidly (waterfall effect) as it catches up with 
the bare dispersion (red solid line) around T. In fact, 
near T, the dressed spectral peak lies slightly below the 
bare band. The width of the spectral function is largest 
in the intermediate energy range of 200-600 meV, where 
its slope also is the largest. This behavior of the spec- 
tral function results from the presence of peaks in the 
real and imaginary parts of the self-energy in the 200- 
500 meV energy range discussed in connections with Figs. 
1(a) and (b) above. It is also in accord with the water- 
fall effect observed in ARPES experiments, although the 
sharpness of the theoretically predicted waterfall in Fig. 
1(c) is less severe than in experiments, which may be due 
to limitations of our model, including the approximations 
underlying our treatment of the susceptibility. 

Fig. 2 considers the case of electron doped (overdoped) 
NCCO. The peak in X' in Fig. 2(a) lies at binding ener- 
gies of 0.5-0.6 eV (in different ^-directions) with a height 
of 0.55-0.7 eV. Correspondingly, the peak in X" in Fig. 
2(b) lies at a binding energy of 0.7-1.1 eV with a height 
of 1-1.4 eV. Comparing these with the results of Fig. 1, 
we see that the self-energy effects in NCCO are much 
larger than in Bi2201. Our computed shift of ^300 meV 
in the position of the peak in X' to higher binding en- 
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FIG. 1: (Color online) (a) and (b): Real and imaginary 
parts £' and S", respectively, of the paramagnon self-energy 
in overdoped Bi2201. Theoretical results at various momenta 
are shown by lines of different colors. The values of £' have 
been shifted by a constant, to produce a zero average value 
of the theoretical E' at the Fermi level. Experimental points 
are from the nodal point for: Bi2201 (gold squares.Ref. (3); 
LSCO (circles, Ref.@); and Bi2212 (diamonds, Ref.|3). Thin 
lines joining experimental points are guides to the eye. Arrow 
marks the location of peak in \" at (f, 0). (c): Spectral den- 
sity in the energy-momentum plane obtained from the dressed 
Green function is shown in a color plot along with the bare 
(red line) and the renormalized (dashed orange line) disper- 
sions. Dots mark the peak positions of the MDC plots of the 
dressed spectral density. 



ergy in going from Bi2201 to NCCO is in good accord 
with the experimentally reported shift of ~300 meV[(|. 
The dispersion underlying the dressed Green function, 
which may be tracked through the yellow dots, is highly 
anomalous and presents a kink-like feature quite reminis- 
cent of the more familiar low ener gy k inks in the 50 meV 
range around the (71", 0)-direction|21|, which have been 
discussed frequently in the cuprates. This strong bosonic 
coupling is also reflected in the fact that the band bot- 
tom in NCCO lies several hundred meVs below the bare 
LDA band in Fig. 2(c). It is interesting to note that 
the self-energies of Figs. 1 and 2 display a 'mirror-like' 
symmetry: The peaks below the Fermi energy in X' and 
X" for Bi2201 in Fig. 1 are smaller than those above the 
Fermi energy, but the situation reverses itself for NCCO 
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FIG. 2: (Color online) Same as the caption to Fig. 1, except 
that this figure shows only the computed results for overdoped 
NCCO. Experimental self-energies are not available and are 
therefore not shown (see text). Solid (dashed) arrow in (a) 
shows location of peak in \" at (^,0) ((it/2, 7r/2)) . 



in Fig. 2 in that now the peaks below the Fermi energy 
become larger than those above the Fermi energy. 

The aforementioned shift of the peak in £' to higher 
energy in NCCO can be understood in terms of the char- 
acteristics of the magnetic susceptibility. Figure [3] com- 
pares in Bi2201 and NCCO the imaginary part which 
is seen from Eq. 2] to be related directly to the real as 
well as imaginary part of the self-energy, x" is seen to 
be quite similar in shape along the T to (ir, 0) line in 
Bi2201 and NCCO, except that in NCCO the band of 
high intensity (the yellowish trace) extends to a signifi- 
cantly higher energy scale. In contrast, x" m the two sys- 
tems differs sharply around (7r,7r). These differences re- 
flect those in the low-energy magnetic response of the two 
cuprates. NCCO with strong magnetic response around 
{■k, 7r) exhibits a nearly commensurate AFM order, while 
Bi2201 is very incommensurate, with peaks shifted to- 
ward (tt,Q). In fact, the high energy peaks in the self- 
energy in Figs. 1 and 2 are tied to the flat-tops near 
(it, 0) at u>i ~ 0.36 eV in Bi2201 (solid arrow), and near 
both (7r,0) at 0J2 — 0.62 eV (solid arrow) and (ir/2,ir/2) 
at w 3 = 0.9 eV (dashed arrow) in NCCO. Above these 
energies the weight in x" falls rapidly, going to zero near 
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FIG. 3: (Color online) Map of the imaginary part of the mag- 
netic susceptibility for (a) hole-doped Bi2201 and (b) electron- 
doped NCCO. Spectral weights are in units of eV~ x . Arrows 
mark the positions of the high spectral weights discussed in 
the text. 



an energy 8t. 

A reference to Figs. 1(a) and 2(a), where the energy 
lui in Bi2201, and the energies u>2 and U13 in NCCO are 
marked by arrows, indicates that the peaks in X' are cor- 
related with these features in the magnetic susceptibility. 
In this spirit, the shift in the peak in £' to higher energy 
in going from Bi2201 to NCCO reflects the fact that fea- 
ture in x" a t ( 7r /2,7r/2) in NCCO (dashed arrow in 
Fig. 2(a)) lies at a higher energy than the (it, 0) feature u>i 
in Bi2201 (arrow in Fig. 1(a)). Notably, when the Stoner 
factor S = 1/(1 — Uxo) is large, a peak in x" arises from 
a peak in Xg( w )i which in turn is associated with nesting 
of features separated by lj in energy. In the present case, 
the nesting is from unoccupied states near the Van Hove 
singularity (VHS) at (it, 0) to the vicinity of the band 
bottom at T, so u)\ ~ 2(i + 2t') - 0.32 eV in Bi2201. The 
larger value of 0J2 in NCCO reflects the shift of the Fermi 
energy to higher energies in an electron-doped material. 

A notable difference between electron and hole dop- 
ing is the low-o; behavior of which is quadratic in 
lu for electron-doping in Fig. 2(b), but nearly linear for 
hole-doping in Fig. 1(b). The linearity for hole-doping, 
reminiscent of marginal Fermi liquid physics, is associ- 
ated here with the proximity of the chemical potential to 
the VHS. This point is considered further in Fig. 4 where 
E" is shown in Bi2201 at the (it, 0) point for three differ- 
ent values of the chemical potential. When the chemical 
potential lies at the VHS (red line), S" varies linearly, 
but when it is shifted by 75 meV above or below the 
VHS, the behavior changes rapidly to become parabolic. 

The strong magnetic scattering discussed in this study 
in the case of overdoped cuprates should persist into the 
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tron and hole-doped cuprates, providing an explanation 
of such effects observed recently in ARPES. Although 
our analysis is limited to the overdoped regime, we ex- 
pect strong magnetic scattering to persist into the under- 
doped regime. This point however bears further study. 

This work is supported by the US Department of En- 
ergy contract DE-AC03-76SF00098 and benefited from 
the allocation of supercomputer time at NERSC and 
Northeastern University's Advanced Scientific Computa- 
tion Center (ASCC). 



FIG. 4: (Color online) Low-energy behavior of £"(a>) at (n, 0) 
in Bi2201 for three different values of the chemical potential 
(j, in relation to the position of the energy Evhs of the van 
Hove singularity. Evhs — ^=0 (red solid line); Evhs — fi= 
+75 meV (green dotted line); and Evhs — n=-75 meV (blue 
dashed line). 



underdoped regime, where the Stoner factor is expected 
to become larger. In fact, this scattering is a precursor 
to the magnetically ordered state near half-filling and 
it is responsible for opening the magnetic gap. In con- 
trast, a number of authors have related the presence of 
waterfall-like effects near half-filling to 'Mott' physics as- 



23|, |24j, but 



sociated with (tt, tt) AFM fluctuations [2z 
have difficulty explaining why these effects persist into 
the overdoped regime. 

The possible doping dependence of U has been an im- 
portant issue in connection with electron-doped cuprates. 
A doping-dependent U is suggested by a number of stud- 
ies in the hole-doped cuprates as well. These include: Op- 
tical evidence of Mott gap decrease [IBj]; ARPES observa- 
tion of very LDA-like bands in optimally and overdoped 
materials; models of the magnetic resonance peak[26|: 
and, a stron gly d oping-dependent gap derived from Hall 
effect studies 27J . The Z 2 -renormalization of %o in Eq. [5] 
bears on this question and gives insight into how the 
value of U enters into the magnetic response of the sys- 
tem. Recall that the susceptibility is often evaluated in 
the literature via Eq. [3] using experimental band param- 
eters, but without the Z 2 factor of Eq. [2] in xo, which 
yields a x scaling ~ Z~ x rather than the correct scaling 
of x ~ This can be corrected by replacing the U in 
the Stoner factor by 



U. 



eff 



= z z u. 



(5) 



Indeed, our Hubbard parameter for NCCO of U = 4i 
is closer to the value at half-filling than is generally 
found. [H 

In conclusion, we have shown that the higher energy 
magnetic susceptibility in the cuprates has considerable 
weight near (it, 0) and that this leads to a high energy 
kink or waterfall-like effect in dispersion in both clcc- 
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